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OVERVIEW 

A  particle  simulation  study  demonstrates  that  an  initial-value  analysis  is 
essential  in  understanding  the  early-time  growth  of  diocotron  modes.  It  is  also 
demonstrated  that  the  equilibrium  properties  of  a  rigid  rotor  (including  rotation 
on  both  the  high  and  the  low  frequency  branches)  are  well  described  by  particle 
simulation. 

NONNEUTRAL  PARTICLE  SIMULATION  CODE 

An  electrostatic,  particle-in-cell  code^*^  is  used  to  model  the  behavior  of  a 
strongly  magnetized,  nonneutral  system  in  the  plane  {x,y)  perpendicular  to  the 
confining  magnetic  field  (along  z).  The  system  is  bounded  by  two  planar  walls 
(located  at  x  =  ±a)  on  which  the  value  of  the  potential  is  zero  (ground).  The 
code  is  2-1/2  dimensional,  i.e.,  it  follows  3  velocity  components  (ua;,Uj„Uj)  and 
two  spatial  coordinates  (x,y)  on  a  time  scale  smaller  than  a  gyroperiod.  Time 
is  scaled  to  the  spatially  averaged  <  The  y  coordinate  is  periodic 

(with  length  L)  and  the  x-variation  is  ignorable.  In  this  study  we  consider  two 
different  initial  configurations:  a  nonneutral  plasma  slab  parallel  to  the  walls, 
and  a  small  rotating  cylinder  far  from  the  walls. 

EVOLUTION  OF  A  NONNEUTRAL  PLASMA  SLAB 

A  pure  electron  plasma  slab  having  uniform  density  and  width  2rp  is  centered 
between  the  walls  at  <  =  0.  The  slab  is  initieilized  by  giving  each  electron  an 
initial  velocity  Vy(i=0)  =  u>^xfflei  corresponding  to  the  expected  x  E.  drift. 
A  small  random  velocity  is  added  to  yield  a  low  temperature  Maxwellian  in  the 
drift  frame.  The  system  is  divided  into  32  grids  along  x  «ind  1024  grids  along  y. 
More  than  32,000  particles  are  typically  followed,  and  the  system  is  well  below 
the  Brillouin  limit;  i.e.,  Qe/<^pe  =  6- 

Figure  la  shows  the  configuration  space  during  the  early  stage  of  the  dioco¬ 
tron  instability  and  Figure  lb  exhibits  the  formation  of  vortices  at  a  later  time. 
The  time  evolution  of  each  Fourier  mode  (having  k„  =  2irnlL)  is  monitored 
using  the  potential  along  x  =  0.  It  is  found  that  there  is  an  inherent  mixture  of 
transient  and  beat  behavior  which  can  be  confused  with  an  instability  or  a  non¬ 
linear  mode-coupling  process.  This  phenomenon  is  illustrated  in  Fig.  2,  where 
the  solid  curves  correspond  to  the  simulation  results.  Figure  2a  shows  that  mode 


20  exhibits  a  monotonic  exponential  growth,  while  Fig.  2b  shows  that  mode  26 
also  grows  initi.  it  is  followed  by  an  abrupt  drop  {t  >  15)  in  amplitude. 

The  initial  grow  Fig.  2b  arises  both  because  of  transients  and  because  of 

the  beat  between  iwo  surface  modes  (one  on  the  top  interface  and  the  other  on 
the  bottom)  moving  past  each  other.  A  proper  description  of  the  phenomenon 
requires  an  initial-value  analysis  too  lengthy  to  report  here.  The  method  uses  a 
Laplace  transform  to  generate  a  differential  equation  in  x  which  is  solved  via  a 
Green’s  function  using  the  method  of  Case^.  The  inverse  transform  is  obtained 
by  evaluating  the  residues,  but  one  remaining  integration  over  x  must  be  done 
numerically.  The  result  of  this  analysis  is  indicated  by  the  dashed  lines  in  Fig.  2. 
It  is  evident  that  mode  20  corresponds  to  a  true  iiistability,  while  the  unusual 
evolution  of  mode  26  is  a  transient  fei.ture  not  contained  in  the  steady-state 
stability  analysis. 

The  growth  rates  of  the  unstabl  ■  observed  in  the  simulation  are  shown 
in  Fig.  3.  These  growth  rates  are  r:.  .  ^d  by  recording  the  potential  energy 

of  each  simulation  mode  over  a  finite  tune  ii  terval.  The  solid  curve  is  the  well- 
known'*  expression  for  the  growth  rate  7  of  an  unstable  mode  (of  wave  number  k) 
of  a  slab  of  width  2rp  in  the  limit  a  — ♦  00; 


Excellent  agreement  between  the  asymptotic  theory  and  the  simulauon  results 
is  obtained  over  a  wide  range.  The  possible  discrepancies  associated  with  the 
high-/?  modes  (18-21)  are  likely  to  be  due  to  the  sampling  of  the  growth  rates 
at  eaxly  times.  The  initial-value  theory  shows  that  early  in  time  these  growth 
rates  are  larger  than  their  asymptotic  values.  We  have  also  simulated  systems 
with  different  ratios  of  rp/a  and  found  excellent  agreement  with  the  theoretical 
predictions.  As  Vp/a  —*  1,  the  system  remains  stable,  as  expected. 

RIGID  ROTOR  EQUILIBRIUM 

We  have  investigated  the  behavior  of  a  small  rigid  rotor  of  radius  tq  bounded 
by  planar  walls  (ro/a  0.16)  as  sketched  in  Fig.  4.  The  magnetic  field  strength 
Bo  has  been  varied  over  the  range  1.5  <  Slg/u;pe  <  3  in  order  to  measure  the 
equilibrium  rotation  frequencies  and  the  radial  oscillation  frequencies.  For  each 
choice  of  Bq  the  system  is  initialized  by  assigning  to  each  particle  an  azimuthal 
velocity  ra>,e,  to  which  a  small  random  component  is  added  to  represent  a  small 
thermal  spread,  and  where 


corresponds  to  the  two  possible  rotation  frequencies  for  a  rigid  rotor  in  vacuum. 
The  average  radius  and  average  rotation  frequency  of  the  cylinder  are  monitored 


as  a  function  of  time.  It  is  found  that  the  rotor  exhibits  long-term  stable  oscil¬ 
lations  in  both  the  mean  radius  and  the  mean  rotation  frequency.  The  top  of 
Fig.  4  shows  that  the  simulated  system  follows  the  classical  rotation  diagram, 
while  the  bottom  shows  that  the  frequencies  of  the  radial  oscillations  follow  the 
solid  curve  predicted  by  the  magnetron  frequency 

^mag  _  ,  _  ^pe 
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